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ABSTRACT 

Context. We have carried out a semi- automated search for planetary nebulae (PNe) in the INT Photometric H- Alpha Survey (IPHAS) 
catalogue. We present the PN search and the list of selected candidates. We cross correlate the selected candidates with a number of 
existing infrared galactic surveys in order to gain further insight into the nature of the candidates. Spectroscopy of a subset of objects 
is used to estimate the number of PNe present in the entire candidate list. 

Aims. The overall aim of the IPHAS PN project is to carry out a deep census of PNe in the northern Galactic plane, an area where PN 
detections are clearly lacking. 

Methods. The PN search is carried out on the IPHAS photometric catalogues. The candidate selection is based on the IPHAS and 
2MASS/UKIDSS colours of the objects and the final candidate selection is made visually. 

Results. From the original list of ~ 600 million IPHAS detections we have selected a total of 1005 objects. Of these, 224 are known 
objects, leaving us with 781 PN candidates. Based on the initial follow-up spectroscopy, we expect the list to include very young and 
proto-PNe in addition to genuine, normal PNe (~16%) and emission line objects other than PNe. We present additional criteria to 
select the most probable PN candidates from our candidate list. 

Key words. Surveys - ISM: planetary nebulae: general - Stars: binaries: symbiotic 

1, Introduction IPHAS is the northern counterpart of the AAO-UKST Ua 

. t ^T^^T X , t . t t n Survey which mapped the Southern Galacti c plane (\ b \ < 10 

Planetary nebulae (PNe) are the evolutionary end products of degrees) and di scovered ~ 1200 new PNe (iParker et aUbooel: 

most low and intermediate mass stars (approximately 1 to 8 ^Miszalski et alj|200i). Both of the surveys are of similar depth 

solar masses). About 2700 PNe have been detected so fom (2-5 Rayleighs) for diffuse, extended emission. However, the 

di^Galaxy, including the objects listed by j Acker et al. | ([19^ resolution and photometric quality of IPHAS is better, making 

T^.^^ TJ^cf ^ fo ^d 1" the Southern Hemisphere by ■ ^ p^^ 

the AAO-UKST Ha Survey jParker et alJl2005t iMiszalski etaP ^ 
l2008h . However, the expected Galactic PN population is much 

larger, with estimates vary ing from 28000 ± 5000 based on I™AS mapped 1800 degrees^ of the Northern Galactic 

observations (iFrew & Parkeii i2006) to 46000 ± 13000 PNe (of ^^^^^ ^^^^ between b= -5 to +5 degrees) m three filters usmg 

size < 0.9 PC) deriv ed from stellar population synthesis models ^^e INT Wide Field Camera at the Observatorio del Roque de los 

(lMoe& De Marcol[2006h . In particular, there is a clear lack of Muchachos (La Palma, Spam). A narrow-band Ua filter (central 

PN detections in the centre o f the Galactic Plane (see for exam- wavelength and width: 6568 A/95 A) and two Sloan filters (/ 

pie Fig. 6c in MiszalskiitiD|2008) due to difficulties in observ- and /') were used for matched 120, 30, and 10 s exposures, re- 

ing nebular objects in areas of high extinction and where confu- spectively. For point sources the survey covers the magnitude 

sion with other types of nebulae (e.g. Hii regions) is significant, range from / ~ 13 to 22 mag (5cr detection limit, magnitudes 

are in the Vega system) with the median magnitude limit being 

Send offprint requests to: K. Viironen /=21 mag. Each IPHAS field is observed twice at two closely 
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overlapp ed pointing s . More information about the survey can be 
found in 'Drew et^ (l2005h . 

To search for new PNe in IPHAS, we use two search meth- 
ods. A semi-automated method was used based on IPHAS pho- 
tometry to find mainly small angular diameter PNe (typically 
< S")? and visual inspection of IPHAS mosaics to find mainly 
extended PNe (typically > 5'')- The two methods overlap in the 
size of objects found, leading to a complete search for PNe with 
diameters from subarcseconds to many arcminutes. Here we 
present the search method and the results of the semi-automated 
PN search. The resul ts of the visual PN s earch will be pubH shed 
in the near future dCorradi et al. 1 120091: ISabin et al.1 120091 S09 
from now on). 

A search for barely resolved or p oint-like PNe based on a 
catalogue o f 4853 Ha em itting stars (^ Witham et al.ll2008h was 
published in! Viironen et a l. ( 2009) (V09 from now on). Here we 
carry out the search in the complete IPHAS photometric cata- 
logue, covering almost the whole Northern Plane including also 
non-stellar objects and objects without detection in the i' band. 
We aim at picking up all the possible PNe included in the photo- 
metric catalogue but due to t he source detection tech nique used 
in the catalogue generation dOonzalez-Solares et al.|[2008h . we 
are more likely to discover compact PNe. In addition to the over- 
all aim of the IPHAS PN project, gaining better understanding 
of the total Galactic PN population, we are especially interested 
in the compact PNe as they are probably either very young or 
distant objects. Thanks to the good spatial resolution of IPHAS, 
we are able to resolve subarcsecond objects which are ideal can- 
didates of young or proto-PNe. These, in turn, give us valuable 
information about the poorly known late stages of stellar evolu- 
tion of the PN progenitor stars. 

The discovery of even a few new very young PNe (vyPNe) 
would be important as only a few of them are known and not 
many are expected in the Galaxy. For example, assuming that 
the PN lifetime is 20000 years and the estimated Galactic popu- 
lation 30000, we can expect that the Galactic population of PNe 
younger than 1000 years is 1500. Of the expected PNe, only 
-2700 are detected and of these about 200 are in the IPHAS area. 
If the constraints of the typical searches for large PNe and for our 
search of vyPNe were the same, we could expect to find only ~ 
10 new vyPNe in the IPHAS area. But as our search is much 
deeper and has better spatial resolution than previous searches, 
we expect to find a larger number of vyPNe. 



2. Search method 

The source of the semi- automated PN search is the IPHAS 
object catalogue, including all the IPHAS observations until 
January 2008 (inclusive), thus covering 93% of the total IPHAS 
area. The initial data release has been published and a detailed 
description of the data reduc tion and products is provided in 
iGonzdlez-Solares et al.1 (l2008l) . 

Our PN candidate search consists of cleaning the catalogue 
data from possible false detections and of choosing the PN can- 
didates due to their loc ation in the IPHAS a nd Two Micron All 
Sky Survey (2MASS, ISkrutskie et"aDl2006h colour-colour dia- 
grams. 

2.1. The catalogue search 

The IPHAS catalogue includes about 600 million detections 
of Northern plane objects, most of them detected two or more 
times. We require detections at least in the Ha and / filters. We 



can assume that a PN is visible in the continuum / filter due to 
presence of the Ha line while a detection in the i' filter is not 
assured. We also require a matching distance between the detec- 
tions in diff'erent filters of < 1.4'' for nebular objects and < V 
for stellar objects. These limits were defined empirically from 
IPHAS detections of known PNe. Based on the known IPHAS 
colours of PNe (see V09) we start our candidate selection by 
making a generous colour cut: / - Ha > 0.25(/ - /') -\- 0.55. 
This cut line was chosen to include all the known PNe in the 
IPHAS two-colour diagram while cutting out the main sequence 
of stars. T he line was aligned along the approximate reddening 
vector (see lCorradi et al.ll200 8'). To remove possible artifacts, we 
discard detections at the borders of the CCDs and in areas of bad 
pixels. In the IPHAS catalogue, all detections are classified as ei- 
ther Saturated, Star, Probable st ar. Probable ex tended. Extended 
or Noise (see Table 7 in Gonza lez-Solares et al. 2008). We re- 
move detections classified as Noise or Saturated, except those 
that are saturated only in Ha in order not to eliminate possible 
very bright Ha emitters. After these initial selection steps, we 
reduced the number of possible candidates to ~ 14 million. 

Most of the possible candidates after the initial selection 
clearly have to be false detections. Studying the initial selection, 
we found the following problems: in cloudy nights, thicker cloud 
cover in the r' than in the Ha image can cause apparent Ha emit- 
ters; areas of diff'use Ha emission contain detections of arbitrary 
nebular objects; in some cases a satellite crossing the Ha-image 
is matched with a star at the same position in / and i' causing it 
to be detected as an Ha emitter. To remove these kinds of prob- 
lems, we require that the objects must be detected at least twice 
and at least one of these two observations must be in a field that 
fulfils the IPHA S quality criteria, in line with th e PhotoObjBest 
requirements in iGonzalez-Solares et al.l (l2008h . About 99% of 
the IPHAS objects are observed at least twice, which allows us 
to apply this selection. For bad quality fields, their observation 
will be repeated later so that the objects possibly lost in this step 
will eventually be recovered. After these selection steps we are 
left with ~ 1.3 million detections. 

So far we have been working with individual IPHAS detec- 
tions. As a next step, we average the magnitudes of the detec- 
tions for each object. For this we use only the magnitudes in the 
fields fulfilling the IPHAS quality criteria. In the case that an ob- 
ject has an /'-measurement in some of the observations but not 
in others, only the good quality images with detection are av- 
eraged. We also make a cut in Ha magnitude, m(Ha)< 19. Our 
motivation for this magnitu de cut is that the limiting magn itude 
of IPHAS is m(/) ^ 20-22 (iGonzalez-Solares et al.ll2Q^ . The 
theoretical maximum IPHAS / - Ha colour of a pure Ha emit- 
ter (assuming that all the flux in the r' filter comes from the Ha 
line) is 3.1 (Drew et al. 2005). Therefore to avoid being biased 
towards objects with smaller Ha excess at faint magnitudes, the 
magnitude cut is a reasonable choice. This also reduces problems 
caused by rising photometric errors at fainter magnitudes. After 
the averaging and the magnitude cut, we are left with ~ 165000 
objects. 

At this stage, the most frequent remaining false detections 
are close optical binaries. Often these objects are not resolved in 
Ha because of the slightly broader PSF in this filter due to the 
longer exposure time. Hence the binary is measured as an ex- 
tended object in the Ha image while it is resolved in / and i\ 
The two stars in / and i' are then matched with the same 'nebu- 
lar' Ha object. We removed these false detections by eliminating 
all objects that have two / detections within T' for the same Ha 
detection. The radius of T' was defined empirically. After re- 
moving this kind of false emitter, we are left with 88000 objects. 
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Fig. 1. Two-colour diagrams - L^/^ IPHAS colours, Right 2MASS colours. Shown are the locations of the IPHAS selected PN 
candidates (diamonds). The ones located also below the cut line in 2MASS are shown as filled diamonds. IPHAS +2M ASS selected 
PN candidates are shown as circles. Main sequence tracks at diff'erent reddenings are shown as solid lines and the colour cuts applied 
in the candidate selection as dashed lines (see text). 



For these objects, the next selection step is based on their IPHAS 
and 2MASS colours. 



2.2. Colour criteria 

After having cleaned the candidates using the above steps, near- 
infrared counterparts of the remaining IPHAS objects were 
sought. A matching radius of 1.4'' was used when search- 
ing for matches in the 2MASS Extended Source Catalogue, 
while for the Point Source Catalogue a T' search radius was 
used. Note that the external astrometric precision of IPHAS 
with respect to 2MASS is g enerally better than 0.1 arcsec 
(iGonzalez-Solares et al ] 120081) . If the object has both an ex- 
tended and a point source counterpart, the latter was selected. 
In the case that no 2MASS counterpart was found, UKIDSSQ 
data was used in the colour selection instead. Having IPHAS 
and 2MASS/UKIDSS information in hand, the PN candidates 
were selected based on their locations in the IPHAS (/ - Ha) vs. 
(/ - /') and 2MASS (/ - H) vs. (H - Ks) diagrams. No colour 
correction was applied to UKIDSS data as the c orrections rela- 
tive to the 2MASS photometric system are small (iHodgkin et al.l 
[2009h . The colour selection is based on the location of known 
PNe and other classes of emission line objects and normal stars 
in thes e diagrams. The diagrams were discussed in Corradi et al. 
(l2008h and V09 and the candidate selection method is the same 
here as in V09. The zone 1 and 2 (see Fig.[T]) introduced in V09 
are defined by the following equations a) and b), respectively: 

a) / - Ha^ > 0.25(/ - i') -h 1.9 or 

r - Ha > 1.9 (for objects without i' magnitude) 

b) 0.25(/ - -h 0.87 < / - Ha^ < 0.25(/ - f) -h 1.9 or 
0.87 < / - Ha < 1.9 (for objects without i' magnitude) 



1 The UKIDSS project is defined in lLawrence et al.l (|2007l). UKIDSS 
uses the UKIRT Wide Field Camera (WFCAM,[C asali et al. 2007) and 



a photometric system described in 'Hewett et al. (200 6). The pip eline 
processing and scie nce archive are described in Irwin et al.l (l2009h and 
iHamblv et al.l (l2008i) . We used data from the 4th data release. 



and 



J-H < IM(H- Ks)-035. 

The objects that fulfil the IPHAS criterion a) and the 2MASS 
criterion b) have the most probable colours to be genuine PNe. 
To be more complete in these highest probability zones, we also 
selected the objects having only a single IPHAS detection for 
visual inspection if they fulfil these colour criteria. 

After the colour selection we are left with 4740 potential PN 
candidates. 



2.3. Visual inspection 

As a final selection step, all 4740 objects were checked visu- 
ally in their IPHAS preview images. With a few exceptions, the 
remaining visually distinguishable mimics are detections in ar- 
eas of large Ha emitting nebulae such as H ii regions or super- 
nova remnants. In some cases, a faint star projected on an ex- 
tended nebulosity can cause a false Ha emitter detection because 
of poor subtraction of the highly variable background. If this is 
clearly the case, these objects are removed, otherwise the object 
was retained. In general, if the non-PN nature was clear from the 
morphology, the object was discarded from the candidate list in 
this visual study. If not, the object was included in the candidate 
list and the final decision on its nature left for future spectro- 
scopic studies. Very extended objects are discarded, as the mor- 
phology is not well seen in the preview images. These objects 
will be recovered in the extended nebula search (S09) of IPHAS 
image mosaics. Generally this size limit is ~r but a few more 
extended nebulae with bright cores are included in our candidate 
fist. 



2.4. Results 

The final result of our PN search is a list of 1005 objects. These 
are listed in Table 1 which is available only electronically. We 
searched first for possible SIMBAD entries for our candidates 
and found 311 matches. Of these, 168 have a primary object type 
of PN or possible PN. Other relatively frequent objects worth 
mentioning are 34 emission objects or emission line stars, of 
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Ha r' 




Fig. 2. Ha, / and i' IPHAS preview images, from left to right, of 
a few FN candidates. The image sizes are T x T, north is up and 
east towards left. The IPHAS detections are shown as circles. 
The three images in the bottom row illustrate the selection of a 
condensation in a new extended FN candidate. 



which three have as a secondary object type PN and according 
to the literature two of these are indeed PNe, 18 Hii regions, 19 
Herbig-Haro (HH) objects, 13 nebulae of unknown nature, 11 
galaxies and 1 1 young stellar objects (YSOs), one being classi- 
fied as candidate only. Our selection also picked out six known 
novae. We have added the SIMBAD information to our final can- 
didate table. We also checked all references in the literature for 
the 311 objects classified in SIMBAD. In 224 cases we consid- 
ered the object classification well established, and we did not 
study these objects any further. 

The 87 objects for which we considered that the object nature 
is not well known are further studied together with the 694 can- 
didates without any SIMBAD classification. Thus, we are left 
with a to tal of 781 PN candidates. Of the 83 candidates from 
the list o f IWitham et aO (|2008) which were already published in 
V09 (as marked in Table 1), all except 10 objects are selected 
here as well. The mi ssing objects were not selected because the 
IWithametal.l(l2QQ8h selection was made before a recent change 
in the Ha filter calibration had been implemented and the objects 
now fall outside our colour cuts or in the fields not fulfilling the 
IPHAS quality criteria. 

In addition to compact candidates, we also found many ex- 
tended objects due to catalogued condensations in larger neb- 
ulae. In these cases, in the final candidate list, the brightest 
(in Ha) catalogued condensation and the corresponding mag- 
nitudes are listed. In addition, the coordinates of the centre and 
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Fig. 3. The S2N2 diagram of ISabbadin et al.1 (1 19771) showing the 
zones for SN Rs and H ii regions, the 0.85 probability ellipse for 
normal PNe (Riesgo & L6pezl l2006h . and the HH zone (V09). 
Over-plotted is our sample of 27 spectroscopically observed PN 
candidates (diamonds). Three new IPHAS vyPNe studied in V09 
are shown as filled diamonds. 



the largest extension of the nebula are given. A detailed study 
of the extended nebulae in IPHAS will be published later (S09). 
Examples of new IPHAS PN candidates are shown in Fig.O 

3. Completeness and reliability of the catalogue 

The catalogue we base our study on covers 93% of the IPHAS 
area. The PN search in the remaining 7% will be carried out 
when the data are available. Our aim is to discover as many new 
Galactic PNe as possible. However, as the search presented here 
is carried out in the IPHAS photometric catalogue, we can aim 
to be complete only for compact PNe (< 5'')- Of these objects, 
we expect to lose 4% in our candidate selection process: 1% due 
to the requirement of a minimum of two detections and 3% due 
to the applied colour selection. 

We tested this in p ractice by c heckin g if the 63 known 
compact (< 5p PN e (lAcker et al.l Il994l: iParker et al.l l2QQ6l: 
MiszalskiWaD 12008") in the IPHAS area and inside the IPHAS 
brightness limits were selected by our algorithm. Six of these 
were not picked up: 4 because of problems in the IPHAS cata- 
logue generation which sometimes splits a single extended ob- 
ject into several objects, and 2 due to our selection algorithm, in 
line with our prediction above. 

We have carried out our search aiming to be as complete 
as possible rather than to have a clean selection of PNe. Of 
the possible mimics, we expect that the most common are T- 
Tauri stars because they enter our colour selection (see V09 and 
Section 14.2.11) and they can be expected to be numerous objects. 
Possible mimics are also symbiotic stars because these objects 
enter clearly into our colour selection (see V09) and their pre- 
dicted number is large (from 3 x 10^ to 4 x 10^ Galactic systems, 
see Corradi et al. 2008, and references therein). Further potential 
mimics are HH objects linked with young stars, although these 
are discarded in the visual inspection if the morphology allows 
a clear separation from PNe. Also, massive young stellar ob- 
jects (mYSOs) and Hii regions are possible mimics. The latter 
were again usually discarded in the visual inspection, however 
compact H ii regions especially might have remained in the final 
candidate list. 
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4. The nature of the candidates 

4.1. Follow-up spectroscopy 

In order to confirm the PN nature of our candidates, optical 
spectra (as a minimum) are needed. For this purpose we have 
started a follow-up spectroscopy campaign. So far we have spec- 
tra for 69 of the PN candidates using the 2.1m telescope at 
the Observatorio Nacional de San Pedro Martir (SPM), Mexico, 
Isaac Newton Telescope (INT), William Herschel Telescope 
(WHT), and Nordic Optical Telescope (NOT), La Palma, Spain, 
and the Australian National University 2.3-m telescope at Siding 
Spring Observatory (SSO), Australia. These include the 19 spec- 
tra mentioned in V09. In this paper, the results from the spectro- 
scopic study are only used to assess the efficiency of our PN se- 
lection. Detailed study of the individual spectra will be published 
in subsequen t papers. A few of t hes e objects have been alre ady 
pubHshed bv iParker et"aD (l2QQ6h or iMiszalski et al.1 (l2QQ8h (as 
indicated in Table 1). However, as they were independently dis- 
covered from IPHAS, they will be included here in our follow-up 
spectroscopy statistics. 



4.1 .1 . New IPHAS PNe, symbiotic stars and emission line 
stars 

Of the 69 observed objects, 38 show a spectrum with nebular 
lines typical of normal PNe. We refer to these as IPHAS nebulae 
in the following discussion. Closer inspection of the 38 spectra 
led us to classify 4 of them as likely symbiotic stars (for de- 
tails a bout recognising symbiotic star spectra, see ICorradi et al.l 
I2QQ8I) . We refer to these as IPHAS symbiotics. 

To analyse the nature of the remaining objects we placed 
them into the log(HQf/[Nii]) vs. log(HQf/[Sii]) (S2N2) diagnos- 
tic diagram (as discussed by VQ9) , which is a revised version 
of the original S2N2 diagram bv .Sabbadin et al.1 (Il977h . 27 of 
the remaining 34 objects show all the necessary emission lines 
in their spectrum, see Fig.[3l Of these, 16 are located within the 
area of the 0.85 PN probability ellipse and outside the overlap- 
ping areas of other types of objects. We will refer to these 16 
objects as IPHAS PNe. 

It is important to emphasise that objects not included within 
the 85% probability ellipse but showing typical PNe lines may 
also very well be genuine PNe. The same can be true for the 7 
IPHAS nebulae not plotted in Fig.[3l As a matter of fact, we have 
marked in Fig.[3]three new PNe discussed in V09 and they all are 
located outside the PN ellipse (one of the V09 new young PNe 
does not show [S ii] lines in its spectrum and is thus not plotted 
here). The objects outside the ellipse are probably peculiar and 
a detailed study is needed in order to classify them. 

The 3 1 clear non-PNe belong to other classes of Ha emitting 
objects such as YSOs. They are interesting newly discovered ob- 
j ects on their own an d a study of their spectra is in progress 
(IValentini et al.ll2QQ9h but in our search for new PNe they are 
set aside as mimics. We will refer to them as IPHAS emission 
line stars, representing the sample of undesired objects in our 
candidate selection. 



4.1 .2. PN discovery rate 

The sample of objects with follow-up spectroscopy is small for 
robust statistics but we estimate that about 16% of the 781 can- 
didates are genuine, normal PNe. We derived this result with the 
following reasoning. 28% of the candidates are located in zone 1 
in the IPHAS two-colour diagram. There are 16 IPHAS PNe, 10 



IPHAS emission line stars and 3 IPHAS likely symbiotic stars 
in zone 1 . From these data we can deduce that 55% of the zone 1 
objects are PNe. On the other hand, from the location of known 
PNe in the IPHAS two-colour diagram (see V09), we know that 
the zone 1 objects represent 95% of the PN population - the rest 
being in zone 2. This leads to the quoted total PN match rate of 
0.28x0.55/0.95 = 0.16. 

The statistics can be biased by the fact that 38 out of 69 
(55%) of the spectra are of objects classified as stellar in IPHAS 
while 582 out of 781 new objects (75%) are stellar according 
to our diameter study (see below). Also, the statistics are based, 
on the one hand, on properties of known PNe and on the other 
hand on properties of new normal PNe (i.e. the group of IPHAS 
PNe, new objects which are located inside the 85% probability 
ellipse of PNe in the S2N2 diagram. Fig. [3]). In addition to this 
kind of normal PNe, we assume that our candidate list includes 
a fair number of young PNe, similar to those presented in V09. 
We note that three out of four of the new young PNe in V09 
are located in zone 2 in the IPHAS two-colour diagram and in 
addition three of them are located outside the 85% probability 
ellipse of PNe, while one of them does not show the [S ii] dou- 
blet in emission in our spectrum. 

However, we note that 90% of our candidates are of diameter 
< 5''. Even a PN match rate of only 16% means that our can- 
didate list includes 112 compact PNe in the Northern Galactic 
plane. This would roughly triple the number of known compact 
PNe in the IPHAS area. Considering that we expect to discover a 
so far hidden population of vyPNe, like the four objects in V09, 
we assume this number to be even higher. 



4.2. Diagnostic diagrams 

In addition to the IPHAS and 2MASS two-colour diagrams, we 
studied the location of our candidates in various diagnostic di- 
agrams combining IR AS (The Infrared Astronomical Satellite, 
iBeichman et aP Il988h 2MASS and MSX (Midcourse Space 
Experiment. lEgan et al.lll999| ) infrared data, in cases where they 
were detected. 



In all of the diagrams we have only plotted those PN can- 
didates which have the needed fluxes well measured. We have 
separately indicated the locations of known PNe included in our 
candidate list. We have also marked the locations of the IPHAS 
PNe, IPHAS emission line stars, and IPHAS symbiotic stars. As 
the sample of new IPHAS objects (for which we have follow- 
up spectroscopy) is small, we have plotted them even when only 
a limit of the flux is available. In these cases, arrows indicate 
the direction for the possible location of the source in the dia- 
gram. As we expect our major mimics to be T-Tauri stars and 
symbiotic stars, we have studied the location of these objects in 
those diagrams where this work has not previously been done by 
other authors. For this p urpose we have searche d the data for the 
symbioti c star sample oflBelczynski et aP (|2000|) and the T-Tauri 
stars of Herbig& Belli (Il988h . 

As mentioned above, in addition to normal PNe, we expect 
that our candidate list includes vyPNe or even proto-PNe, like 
the four objects presented in V09. To get an idea about the ex- 
pected locations of these kinds of objects in the diagnostic dia- 
grams, we have also plotted in them the 'very-likely post-AGB 
stars' of ISzczerba et al.l (l2007h . 
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Fig. 4. L^/^ Location of the new IPHAS PNe (filled diamonds), new IPHAS emission line stars (fi lled circles), and new IPHAS 
symbiotic stars (filled triangles) in the IP HAS two-colour diagra m. Plotted are also the locations of iHerbig & Bdll (flOSSk T-Tauri 
stars (squares) and the post-AGB stars of ISzczerba et alJ (l2007l) (plus signs). The known PNe in our candidate list are shown as 
asterisks. Main sequence tracks are shown as solid lines and the colour cuts applied in the candidate selection as dashed lines (see 
text). Right 2MASS two-colour diagra m. The same symbo ls and line- styles are used as in the left hand panel. T-Tauri stars are not 
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plotted as their location was studied in lCorradi et al.l (l2QQ8h . 



4.2.1. IPHAS and 2MASS 

The IPHAS and 2MASS two-colour diagrams for the objects for 
which we have spectra are shown in Fig.H] The error bars are not 
plotted but are given in Table 1 . However, generally the errors 
are about the same size as the symbols used. The locations of 
the newly discovered objects are in line with the locations of 
the previously known objects (see V09). The location of T-Tauri 
stars in IPHAS was not shown in V09 and here we see that some 
T-Tauri stars enter our zone 2 in the IPHAS diagram. 

Only 42 post-AGB stars fromlSzcz erba et al.l (12007') have an 
IPHAS counterpart. These objects do not necessarily yet show 
Ha in emission which is apparent in the IPHAS two-colour dia- 
gram. Below our selection line in 2MASS, the post-AGB objects 
are located towards higher J-H and H-Ks colours as compared 
to PNe. We note that the new young V09 PNe also generally have 
lower / - Ha colours and higher J-H and H - Ks colours than 
the bulk of the normal PNe. 

To conclude, we consider a higher probability for our candi- 
dates to be normal PNe if they are located high up in the IPHAS 
two-colour diagram and towards low J-H and intermediate 
H - Ks colours in the 2MASS two-colour diagram. However, 
for vyPNe we expect a lower position in the IPHAS two-colour 
diagram and larger J-H and H - Ks colour values. 

4.2.2. IRAS 

We cross-matched our candidate list with the IRAS catalogue, 
searching for possible matches in the IRAS Point Source 
Catalog, IRAS Faint Source Catalog, IRAS Cataloged Galaxies 
and Quasars, IRAS Serendipitous Survey Catalog and in the 
IRAS Small Scale Structure Catalog. When searching for IRAS 
data, a W search radius was used. This choic e is based on the 
pos itional accuracy of extended IRAS objects (iBeichman et al.l 
11988). Of our 781 PN candidates, 51 are well measured (at least) 
in the 12yum, 25yum and 60yum IRAS bands. 

IPottasch et al.l (Il988l) defined diff'erent areas marking the 
location of Hii regions, OH/IR stars, PNe, normal late 



-0.5 



Galaxy 




du|ftQ>W 



%~3 ^ ^ i pn" 



— ■^ :/-vTT/T-n+ i +"^.4-+ ,J^J + ^ W_ ^ + 



- * 



-0.5 



0.5 

log(F25/F12) 



1.5 



Fig. 5. IRAS two-colour diagram. The symbols used are the 
same as in Fig. |4l Shown are also the locations of our candi- 
dates (open diamonds) and the iBelczynski et al.l (2000) symbi- 
otic stars with an IRAS couterpart (open triangles). 



type stars, and galaxies in the IRAS F(12yum)/F(25yum) vs. 
F(25yum)/F(60yum) diagram, as shown in Fig. [5] We also added 
to the diagram the locations of T-Tauri stars as studied by 
Emerson (1987) and of outflow source s (both young and e volved 
stars with fast bipolar outflows) from lZijlstra et al.r(l200l|) . 

The locations of the newly discovered objects are again in 
line with expectations except for one IPHAS PN which is clearly 
above the PN area, in the Galaxy/Hii region zone. However, the 
location of this object, IPHASXJ195248.8-h255360,in the S2N2 
diagram and its morphology (see Fig. [7]) clearly point to a PN 
nature. We consider a higher probability for our candidates to be 
a PN if they are located inside the PN box. But we also point out 
that a location outside of the PN box does not automatically rule 
out a PN nature. 
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4.2.3. MSX + 2MASS infrared diagrams 

We also searched for MSX counterparts for our PN candidates. 
A 5'' search radius was used based on the MSX positional accu- 
racy dEg an et al. 1999). 122 of our PN candidates have an MSX 
co unterpart well ni e asured at least in one of the MSX filters. 

iLumsden et al.l (l2QQ2h studied various diagnostic diagrams 
combining MSX and 2MASS flux ratios in order to select mYSO 
candidates. For this purpose, they placed Herbig Ae/Be stars, 
mYSOs, methanol maser sources with radio emission (which 
trace massive star forming regions), compact H ii regions, carbon 
stars, OH/IR stars, and PNe into the diagrams. We refer here to 
the diagrams shown in Figures 3, 4, 9 and 10 in lLumsden et al.l 
(2002) as LI, L2, L3 and L4, respectively. In each of these dia- 
grams we defined a line along the reddening vector which lim- 
its the zone where most of the known PNe in the diagrams are 
located. These cut lines in the LI, L2, L3 and L4 diagrams, re- 
spectively, are as follows: 

log[F(14yum)/F(8yum)] > -6Mlog[F(l2/im)/ F(U/um)] -h 2 

log[F(Ujim)/F(l2jim)] < 0.561og[F(21yum)/F(8yum)] - 0.21 

log[F(^/um)/F(Ks)] < 3.031og[F(12yum)/F(8yum)] - 0.86 

log[F(Ks)/F(J)] < 29.29 log[F(21yum)/F(12yum)] -4.82. 

We studied the locations of our candidates, symbiotic stars, 
T-Tauri and post-AGB stars in the diagrams LI, L2, L3 and L4 
and found that all of them separate PNe from T-Tauri and sym- 
biotic stars quite nicely. However, we note that in all these dia- 
grams, the PN zone overlaps with the zones occupied by com- 
pact H II regions and methanol maser sources with radio emis- 
sion. There is also a partial overlap with the OH/IR star zone. 
In the L3 and L4 diagrams, also some Herbig Ae/Be stars and 
mYSOs enter the PN zone. Nevertheless, as we expect T-Tauri 
and symbiotic stars to be the most frequent mimics in our cata- 
logue, we consider these diagrams useful in weighting the prob- 
ability that our candidates are a genuine PNe. However, the post- 
AGB stars are not separated from the rest of the objects in any of 
the diagrams. This can mean that the vyPNe would not be either. 

4.2.4. Discussion on the diagnostic diagrams 

We have studied in various diagnostic diagrams the locations of 
our candidates and the expected most frequent mimics, na mely 
T-Tauri stars and symbiotic stars. Thel Lumsden et al.l ("2002) di- 
agrams do not separate PNe from compact H ii regions but the 
separation from other possible mimics is generally good, espe- 
cially in the IRAS diagram and in the LI and L2 diagrams. The 
IR data certainly can be used to further restrict the possibility 
that our candidates are genuine PNe in the cases where the nec- 
essary IRAS, MSX and 2MASS data for the candidate exist. 

We have also placed post-AGB stars in these diagrams as 
we expect that our candidate list includes, in addition to normal 
PNe, vyPNe or proto-PNe like those presented in V09. Based 
on the locations of the post-AGB stars as compared to PNe we 
expect that the vyPNe are found outside the most typical PN 
zones in IPHAS and 2MASS two-colour diagrams. The results 
in V09 support this observation. In these zones the confusion 
with mimics is higher and thus makes the discovery of vyPNe 
more difficult. In the IRAS and MSX diagrams the locations 
of post-AGB stars overlap with both PNe and normal late type 
stars. This might be an evolutionary trend, the most evolved post- 
AGB stars being located in the PN zone while the less evolved 



ones are located in the zone of late type stars. We find that we 
cannot separate the vyPNe from the normal late type stars based 
only on their infrared colours. Lastly, we note that in the IRAS 
diagram, the distribution of our candidates clearly diff'ers from 
the distribution of known PNe and post-AGB stars. Therefore, 
we conclude that most of our candidates with an IRAS coun- 
terpart are rather mimics than young PNe. We have included all 
available IR data in Table 1 . 



5. Additional data 

5.1. Radio 

Candidates with PN-like colours in the diff'erent optical/IR di- 
agrams, showing an adequate morphology and associated with 
radio continuum sources, are almost certainly PNe. We searched 
for matches to our candidates b y comparing the IP HAS and 
NVSS (NRAO VLA Sky Survey, 'Condon et al."'1998) radio im- 
ages. 63 of the 781 IPHAS PN candidates were found to have an 
1 .4 GHz radio counterpart. We add the radio information to our 
candidate table. 



5.2. Diameters of the candidates 

Normal PNe are expected to be more extended than T-Tauri stars, 
symbiotic stars, post-AGB stars and other stars and we mea- 
sured the diameters of all our candidates. For the 75 clearly ex- 
tended candidates, the diameters and coordinates of the centre 
were measured in the IPHAS Ha images using the FITS viewer 
ds50. For the 706 candidates that appear compact, we have mea- 
sured their full width half maximums (FWHMs). For compari- 
son, we measured the FWHMs of field stars in an area of 10' 
surrounding the objects, selecting the stars within ±0.1 magni- 
tudes from our candidate, or as a minimum, 50 stars as close 
in magnitude to our object as possible. A 5% trimmed and 2-cr 
clipped mean of these defines the stellar FWHM and the RMS 
its error, as well as the error in the FWHM measurement of 
the object. We decided to use both trimming and sigma clip- 
ping for a robust measurement of the stellar FWHM. In order 
to claim that one candidate is resolved, we use a conservative 
limit of five sigma, i.e. we consider that the object is resolved if 
FWHM(object) - FWHM(stars) > 5 x err_FWHM(object). 
Of the 706 compact objects, 124 turned out to be resolved us- 
ing this criterion. Of these, 25 have a deconvolved diameter, 
^JFWHM(object)^ - FWHM(stars)^, smaller than 1 arcsec. 

120 of the resolved candidates are of size < 5'', the size range 
for which our search method is optimised. At diameters greater 
than 5'\ only a fraction of the nebulae where the semi- automated 
algorithm detected a small condensation is found. We hardly find 
any candidates with diameters smaller than 0.5'', as expected for 
instrumental and observational reasons. There are 75 candidates 
with diameters in the range 0.5" - 2". This is the diameter range 
where young/compact/very distant PNe are expected, possibly 
accompanied by a small fraction of nebulae around T-Tauri stars 
and symbiotic stars. From 2" on, the number of objects clearly 
drops with increasing diameter. This can be attributed to the dif- 
ficulty in detecting lower surface brightness objects at large dis- 
tances due to interstellar extinction. 



^ http ://hea- w w w. harvard. edu/RD/ds9 
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Fig. 6. IPHAS Ha images of the 14 extended candidates with PN-like morphology and which have an entry in SIMBAD but whose 
nature is not confirmed. 



6. Extended candidates 

Although our search method is optimised to discover compact 
(<5'0 PNe, we also found 81 extended objects with a diame- 
ter >5''. For these clearly extended objects their morphology 
can help in determining their nature. 36 nebulae show regular 
morphology so that from their appearance alone we can assume 
a higher probability that th ese objects are genuin e PNe. Two 
of these were discovered by iMiszalski et al.l (l2Q08l) . We present 
here the images and short descriptions of the remaining 34 ob- 
jects. Of these, 14 are known in SIMBAD but their classification 
is not well established and they deserve further study. 

6.1. Candidates classified in SIMBAD 

IPHASXJ022045.2+631135 Classified as a galaxy ("Hau et al." 
fT995h . However, it is hardly visible in the IPHAS continuum 
filters, which is an argument against the galaxy classification. 
Possibly an elliptical PN. 

IPHASXJ052708.2+383114 Classified as a reflection nebula 
(|Magakianll20Q3i) but no detailed study of the object is available. 
Possibly a round PN. 



IPHASXJ055242.8+262116 A CO source (IWouterloot et al.l 
119931) . Possibly a round PN. 

IPHASXJ190438.7+0 21424 Cata l ogued as an Ha-emission 
star (Kohoute k & Wehmevejll999l) . ICondon et aP ([1999) have 
obtained an optical spectrum of this object and classify it as 
a PN. Preliminarily classified as a symbiotic star based on a 
MASH spectrum but more observations are needed. An elliptical 
PN or symbiotic star. 

IPHASXJ190954.7+120455 Classified as a possible PN 
( Kronberger et al. 2006) but no spectroscopic confirmation is 
available. Possibly a round or bipolar PN with a brightened 
waist. 

IPHASXJ192751.3+1 40127 Classified as a possible PN 
(|Preite-Martinezl Il988|) but no spectroscopic confirmation is 
available. Possibly a round PN. 

IPHASXJ193 141.3 +224339 Classifications of PN 
( Condon et al.l 1 1999[) , Gala xy (iRoman et al.ll2000h and possible 
PN Preite-Martinez but no spectroscopic confirmation 

is available. Hardly visible in the IPHAS continuum filters and 
thus unlikely to be a galaxy. Possibly a round/slightly elliptical 
PN. 
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Fig. 7. IPHAS Hof images of the 20 extended candidates with PN-like morphology. 



IPHASXJ194040.3 +293010 Classifie d as a possible PN by 
iRoman et all (l2QQd) and as a PN by Kro nberger et alJ (l2006h but 
no spectroscopic confirmation is available and hardly visible in 
the IPHAS continuum filters. Possibly a round PN. 
IPHASXJ194226.1+21 4522 Classified as a possible PN 
(iKronberger et al.l r2006h but no spectroscopic confirmation is 
available. Possibly a bipolar/quadrupolar PN. 
IPHASXJ195209.1+271831 Classified as PN (PN K 3-48). 
Many stu dies refer to this ob ject as a PN but is listed as an Hii 
region by I Acker et al.l (Il994l) . No spectra available in the litera- 
ture. Possibly an elliptical PN. 

IPHASXJ200018.7+365935 Classified as a possible PN 
(lAcker et al.1 1 19941) but with no spectroscopic confirmation. 
Possibly a round PN. 



IPHASXJ200457.3+311416 Listed as a 'Bright Nebula' 
( Kronberger et al. 2006) . Possibly a round or a bipolar PN with 
brightened waist. 



IPHASXJ202549.3+393020 Catalogued as an Ua emission fine 
star ( Kohoutek & Wehmever 1997), PN (Manchado et al. 1989) 
and H ii region ( Garcia-La rio et al.l 1 1997) but with no spectro- 
scopic confirmation available. Possibly a bipolar PN. 



IPHASXJ232713.1+650923 Classified as possible PN 
(iKronberger et al.ll2006 ) but with no spectroscopic confirmation 
available. Possibly a round PN. 
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6.2. Candidates without any SIMBAD classification 

IPHASXJ002059.6+623937 A bipolar nebula. 
IPHASXJ025658.2+585548 Possibly bipolar but a deeper im- 
age is needed to define the morphology. 

IPHASXJ185309.4+075241 A roundish/sHghtly elliptical neb- 
ula. 

IPHASXJ191345.6+174752 A round nebula. Brighter rims, es- 
pecially the western one. 

IPHASXJ191445.1+133220 Round nebula. Brightened rim. 
IPHASXJ191727.0+082036 Faint nebula. Looks round but a 
deeper image is needed. 

IPHASXJ192140.4+155354 Butterfly shaped nebula. 
Preliminarily classified as a symbiotic star based on a 
MASH spectrum but more spectra are needed to confirm the 
classification. 

IPHASXJ192902.5+244646 A round nebula. 
IPHASXJ193532.1+112116 A round nebula. 
IPHASXJ193718.8+202102 A bipolar nebula. Only two irreg- 
ularly shaped lobes visible. 

IPHASXJ193949.8+264047 Shape unclear but if the bright 
central part is the waist of a bipolar nebula, the real extent can 
be larger than the W listed in Table 1. 

IPHASXJ194648.2+193608 A bipolar nebula. The two lobes 
are combined with an elliptical ring. 

IPHASXJ194940.9+261521 A bipolar nebula. The waist is 
clearly brighter than the lobes. 

IPHASXJ195126.5+265839 Faint round/elliptical nebula. 
IPHASXJ195248.8+255360 Bipolar nebula. The waist is bright 
while the lobes are very faint. 

IPHASXJ204118.3+471011 The morphology is not very clear 
but looks like a bright waist of a bipolar nebula with faint lobes 
only partially visible. 

IPHASXJ212000.1+514106 An elliptical nebula. The NW and 
SE rims are brightened. 

IPHASXJ212335.3+484718 A bipolar nebula, bright waist and 
only partly visible lobes. 

IPHASXJ212608.4+542015 An 'S' -shaped nebula. 
Quadrupolar? Deeper imaging needed to define the mor- 
phology. 

IPHASXJ214032.5+564752 A round nebula. 



7. The IPHAS PN candidate list 

We have collected all the information discussed above into a list 
of IPHAS PNe candidates (provided in electronic format only). 
The table contains all the 1005 objects picked up by our selection 
algorithm. The columns of the table are as follows: 
Name: The lAU-registered sexagesimal, equatorial position- 
based source name in the form: IPHASJhhmmss.ss-hddmmss.s 
for stellar and IPHASXJhhmmss.s-i-ddmmss for extended ob- 
jects. "J" indicates the position is J2000. 

RA, DEC: The right ascension and declination (J2000) of the 
object catalogued by IPHAS [deg]. 

Gal: The galactic coordinates of the object in the form 111.1-hbb.b 
where 1 and b stand for galactic longitude and latitude, respec- 
tively, truncated to one decimal place [deg] . 
r, r-i, r-Ha: The IPHAS / band magnitude, / - i' (if measured 
in /') and / - Ha colours of the object. In the case of extended 
objects, the IPHAS magnitude/colour information corresponds 
to the brightest catalogued Ha condensation of the object. 
f_r-Ha: IPHAS photometry flag, meaning that the object is satu- 
rated in Ha. 



cRA, cDE: The right ascension and declination (J2000) of the 
measured centre for the object [deg]. This information is only 
presented for extended objects. 

IRCat: The near-IR catalogue (2MASS PSC, 2MASS XSC, or 
UKIDSS) from which the near-IR data of the object was taken. 
J, J-H, H-Ks: / band magnitude, J - H and H - Kg colours of 
the object. 

IRASCat: The IRAS catalogue (PSC/SSC/SSS) from which the 
IRAS data of the object was taken. 

fl2, f25, f60, flOO : IRAS fluxes at 12, 25, 60 and 100 yum [Jy]. 

fMSXa, fMSXc, fMSXd, fMSXe: MSX fluxes at 8.28, 12.13, 

14.65 and 21.34 yum [Jy]. 

fNVSS: Radio flux at 1.4 GHz [mJy]. 

sFWHM: The stellar FWHM [arcsec]. 

oFWHM: The object FWHM [arcsec]. The FWHM informa- 
tion is provided only for the objects appearing quasi- stellar in 
IPHAS. 

Diam: Diameter for the clearly extended candidates [arcsec]. 
PNIRAS: In the PN zone in the IRAS diagram [yes/no]. 
PNLl: In the PN zone in the LI diagram [yes/no]. 
PNL2: In the PN zone in the L2 diagram [yes/no]. 
PNL3: In the PN zone in the L3 diagram [yes/no]. 
PNL4: In the PN zone in the L4 diagram [yes/no], 
score: The PN score as defined below. 

LSIMBAD: Flag indicating that we consider the SIMBAD type 

not well established. 

type: The types Hsted in SIMBAD. 

SIMBAD: The SIMBAD name. 

crossref: Cross reference with the other IPHAS and AAO- 
UKST Ha Survey selected catalog ues: C = IPHAS symbiotic 
star candidates ( Corra di et aLl l2008V. S = IPHAS extended nebu- 
l ae catalogue (S09), M = PN in M ASH o r MASH-II catalogues 
(iParker etaLlt 2006: Miszalski et al.l l2008i) . W = IPHAS emitter 
catalogue ('Witha m et al.. .2008). 

class: Preliminary classification based on follow-up spec- 
troscopy. A detailed analyses of the objects and the spectra will 
be published elsewhere. Objects with a MASH-II classification 
have been followed up by the MASH project. (IPHAS = spec- 
trum taken but not yet studied, IPHAS PN/MASH-II PN = likely 
planetary nebula, IPHAS nebula = possible PN, IPHAS/MASH- 
II star = Ha emission line star, IPHAS/MASH-II symbiotic = 
likely symbiotic star. See S09 = will be published in S09, See 
V09 = was studied in V09). 

In addition the magnitude and colour errors are given in the 
columns e_50 and the flux qualities in the columns q-50 for the x 
data. 

The PN score is defined from the source location in the 
IPHAS, 2MASS and other IR diagnostic diagrams in the fol- 
lowing way: In the IPHAS two-colour diagram, the distance in 
magnitudes to the lowest selection line is calculated. To this we 
add the calculated distance in magnitudes from the selection line 
in the 2MASS two-colour diagram. If the object is located to the 
left of this selection line the distance is simply set to zero. If 
in addition the object is located in the PN zone in any of the 
infrared diagrams, its score is increased by one for each IR dia- 
grams where the object is in the right zone. So, in summary, the 

^ For the near-IR data an error value of zero indicates that the magni- 
tude, or one of the m agnitudes defining the colour, is a 95% confidence 
upper limit dCutri et al. 2003). 

4 For IRAS PSC and SSC: 3=high quality, 2=moderate quality, 
l=upper limit. For IRAS SSS: A=high quality, B=intermediate qual- 
ity, F=low quality ( Beichman et al. 1988). For MSX: 4 = excellent, 3 = 
good, 2 = fair, 1 = limit (.Egan et al...l999l) . 
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higher the PN score of the object is, the higher we consider the 
possibility that this candidate is really a PN. For example, for 
the known PNe in the candidate list, the average PN score is 3.6 
while for the other types of objects classified in SIMBAD the 
average is 2.5. For the new IPHAS PNe the average score is 3.3 
and for the new IPHAS emission line stars 1.9. 

However, the PN score alone cannot obviously be used to 
discard the possibility that a candidate is a PN. For example, the 
PN score of the known, extended PN, PN G036.9-02, is only 
0.52. Its Ha excess is low due to an underlying star which was 
picked up by our selection algorithm and it is not an IRAS nor 
MSX source which explain the low score. 

In addition to the PN score, the radio and diameter data give 
additional constraints on the probable nature of the object. 

8. Summary and discussion 

IPHAS provides a powerful database for PN searches in the 
Northern Galactic Plane. Thanks to its high resolution and depth, 
especially the two extremes of the PN evolution, young, compact 
PNe and old, faint, and extended PNe are explored. The visual 
search, optimal for the discovery of very extended PNe, is un- 
derway (S09) while here we have presented the semi-automated 
search for PNe based on the IPHAS photometric catalogue. This 
method is optimised to discover small- sized PNe but can also 
discover some extended PNe. In addition to normal PNe, we ex- 
pect that our candidate list includes vyPNe and proto-PNe. 

We present a list of 1005 objects selected from the initial 
catalogue of ~ 600 million IPHAS detections. The SIMBAD 
database and literature study reveals that 224 of these are ac- 
tually known objects, 149 being confirmed PNe. This leaves us 
with 781 PN candidates. About 25% of the selected candidates 
are slightly extended and the clearly extended ones (> 5'') over- 
lap with the visual search by S09 for extended PNe. This ensures 
that the combined IPHAS PN search is as complete as possible. 

When selecting the PN candidates, we have aimed to be com- 
plete rather than having a clean selection of candidates. The first 
results from our follow-up spectroscopy reveal that 38 out of 69 
candidates observed show nebular lines, whereas the remaining 
candidates are Ha emission line stars, most probably YSOs. In 
order to help in the selection of good PNe candidates we have 
studied the locations of the candidates in a number of IR diag- 
nostic diagrams. Further, we provide additional information like 
the sizes of the objects and the available radio data for them. 

The probability that a candidate is a genuine, normal PN is 
higher if it is located in the PN area in the IPHAS and IR dia- 
grams and if it is extended in Ha. All this information is pro- 
vided in our candidate list and the colour-colour diagram loca- 
tions are used to define a PN score whose value is higher for 
candidates with properties like the known, normal PNe. 

However we would like to point out that exploring the ob- 
jects not fulfilling the "normal PN" criteria listed above can be 
surprisingly productive in discovering objects in rare PN phases, 
as is shown by the very young PNe discovered and discussed in 
V09. 
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